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ABSTRACT 

By utilizing satellite-based estimations of the distribution of clouds, we have studied Earth's 
large-scale cloudiness behavior according to latitude and surface types (ice, water, vegetation and 
desert). These empirical relationships are used here to reconstruct the possible cloud distribution 
of historical epochs of Earth's history such as the Late Cretaceous (90 Ma ago), the Late Triassic 
(230 Ma ago), the Mississippian (340 Ma ago), and the Late Cambrian (500 Ma ago), when the 
landmass distributions were different from today's. With this information, we have been able to 
simulate the globally-integrated photometric variability of the planet at these epochs. We find 
that our simple model reproduces well the observed cloud distribution and albedo variability 
of the modern Earth. Moreover, the model suggests that the photometric variability of the 
Earth was probably much larger in past epochs. This enhanced photometric variability could 
improve the chances for the difficult determination of the rotational period and the identification 
of continental landmasses for a distant planets. 



Subject headings: Astrobiology — atmospheric effects — Earth — planets and satellites: atmospheres 
— planets and satellites: general — planets and satellites: surfaces 



1. INTRODUCTION 

Since t he discovery of t he fir st extrasolar 
planet by iMavor fc Queloz (1995), more than 
500 planets orbiting stars other than the Sun 
have been detected using a variety of tech- 



that in the near future it is very likely that Earth- 
size planets can be discovered in large numbers 
(|Lindensmithll2Q03l: [Schneider et al.ll2009[) . 

The atmospheric characterization of giant 
(Jupiter-like) pla nets has alre a dy st arted with 



niques ( Beaulieu et al. 2006t Beuzit et al. 2007 



promising results (|Tinetti et al.l (|2010f ). and refer 



Charbonneau et al.|[2007t lUdrv et al.ll2007l ). Most 
of them are gas giants, as larger planets are easier 
to detect than smaller, rocky planets. However, 
tens of planets in the sup er-Earth mass range 



have also bee n disco v ered (IBeaulieu et al.l 2006 : 
Bennett et all bOOSt ICharbonneau et al.l boOS : 
Mavor et al.ll2009t iQueloz et al.ll2009l ). Ahhough 
we are not yet capable of detecting and explor- 
ing planets like our own ( the smallest planet y et 



known has Rp = 0.127i? r: Ifiatalha et al.1 (|201l[ )). 
ambitious ground and space based projects are al- 
ready being planned for the next years, suggesting 
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ences therein) , and several atmospheric models are 
being developed to help us understand these ob- 
servations. In this sense, it is clear that the explo- 
ration of our own solar system and its planets will 
provide a useful opportunity for method valida- 
tion, enabling more accurate determinations, and 
characterization of extrasolar planets. In particu- 
lar, observations of the solar system rocky planets, 
including Earth, will be key for the search for life 
elsewhere. 

In the last years many studies, both obser- 
vational and theoretical, have been performed 
to investigate how the Earth would look like 
to an extrasolar observer. One of the observa- 
tional approaches has been to observe the light 
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reflected by Earth, referred to as Ear thshine, via 



the dark side of the moon at visible (jPalle et al 
[200^ ,Qiu ct al. 2003) and near-IR wavelengths 
(jrurnbull et alj l2006t) . Another approach has 
been the analysis of Earth' s light curve s from 



remote sensing plat forms (jCowan et al.l 12011 



Robinson et al. l201l[ l . Several authors have also 
attempted to model the diurnal photometric 
varia b ility on an Earth-lik e planet JPo rd ct al] 



20011 iTinetti et all l2006allbl: iFuiii et al 



while other authors such as Woolf et al 



Arnold et al.l (120021 120091) . ISeaeer et al 




and iMontafies-Rodrfguez et al.l (|2005l l2006[ ) have 
attempted to measure the characteristics of the re- 
flected spectrum and the enhancement of Earth's 
refle ctance at 700 nm d ue to vegetation. In addi- 
tion, iPalle et al.l (|2008l ) determined that the light 
scattered by the Earth as a function of time con- 
tains sufficient information, even with the presence 
of clouds, to accurately measure Earth's rotation 
period. Even crude reconstruction of the conti- 
nental distribution could be attemp ted given suf 



ficient signal-to-n oise observations (jCowan et al 
20091 : lOaklev fc C ash 2009). 



However, it is unlikely that, even if we were 
to find an Earth-twin, that planet will be at a 
similar stage of evolution as the Earth is today. 
On the contrary, extrasolar planets are expected 
to exhibit a wide range of ages and evolutionary 
stages. Because of that, it is of interest not only to 
use our own planet, as it is today, as an exemplar 
case, but also at different epochs. 

In this paper, we have used information about 
the surface properties and continental distribution 
of the Earth with the aim of studying the behav- 
ior of the large-scale cloud patterns. We have ob- 
tained empirical relationships between the amount 
of cloudiness and their location on Earth's sur- 
face depending on latitude and underlying surface 
types. This relationships have been used to re- 
construct the possible cloud distribution for dif- 
ferent epochs of the Earth such as 90, 230, 340, 
and 500 Ma ago. With this information, we here 
attempt to reconstruct and understand the pho- 
tometric variability of these epochs according to 
their different geographic and clouds distribution. 



2. CLOUD AND CONTINENTAL DIS- 
TRIBUTION DATA 

The International Satellite Cloud Climatol- 
ogy Project (ISCCP) was established in the early 
1980s. One of its main goals has been to obtain 
more information about how clouds alter the ra- 
diation balance of the Earth. To this end ISCCP 
has been collecting and analyzing satellite radi- 
ance measurements to infer the global distribu- 
tion of clouds, their properties, and their diurnal, 
seasonal, and interannual variations. 

In our analysis we have used the monthly mean 
(D2) fractional cloud cover, derived from a com- 
bination of both VIS and IR measurements, for 
the 23 year period 1984-2006. The data are given 
for 2°. 5x2°. 5 grid cells in latitude and longitude. 
The ISCCP datasellll provides information on sev- 
eral surface, atmospheric and cloud parameters at 
each grid point. In particular total cloud amount, 
and several cloud types are given at each grid 
point. For further d etails on the ISCCP data see 
Rossow et all (|l996[ ). 

Detailed information on the varying distribu- 
tion of continents, deserts, vegetation and ice 
in different epochs of Earth's history, is avail- 
able from Ron Blakey's web pag^ on paleocli- 
mate reconstructions. These geological maps have 
been constructed by using interpretations from 
a wealth of geologic literature and publications, 
avai lable on the project 's web page (see, for exam- 
ple, iFrisch et al.l (j201ll) , and references therein) . 

In order to classify the correspondences be- 
tween the land areas of these paleomaps and the 
ISCCP land types classification scheme, w e used 



an in tensity thresholding method ( Son ka et al 
19931 ). Next, the information was regridded to 



a geographical resolution of 144 (longitude) by 72 
(latitude) cells, equal to the ISCCP cloud prod- 
ucts, as input for our albedo models. Figure [T] 
shows the four historical periods that we have 
chosen to carry out our study. These epochs are 
90, 230, 340 and 500 Ma ago. Note that we have 
assumed that the Earth 500 Ma ago was entirely 
desert as the development of advance plants is be- 
lieved to have taken place in the Late Ordovician 
(450 Ma ago) (|Grav et al.lll985l ). ahhough fungi. 



-'http:/ /isccp. giss.nasa.gov 
^http: / /jan. ucc.nau.edu 
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algae, and lichens might have greened many land 
areas some time before. Whether this is really the 
case for the Late Cambrian or not is not the real 
issue, as our supposition is meant to illustrate the 
albedo properties for a planet with a continental 
crust but not life in its surface, which was prob- 
ably the case for Earth in many previous epochs, 
and a plausible scenario for Earth-like extrasolar 
planets. 

The selection of these four geological epochs 
corresponds to a period of time during which plate 
tectonics has radically changed the face of the 
planet. However, albeit with large temporal excur- 
sions. Earth's atmosphere can be considered sim- 
ilar t o the present f)Hart 1978.: Kasting fc SiefertI 
20021 ). On average the same atmospheric composi- 
tion and mean averaged temperature have existed 
during this period. This is a necessary condition 
for our supposition of cloud distribution related 
to latitude and surface type to hold as a valid ap- 
proximation. 

3. A SEMI EMPIRICAL MODEL FOR 
CLOUDS 

When attempting to study the reflectance prop- 
erties of the Earth in the past, one of the major 
difficulties is the complete lack of reliable informa- 
tion on cloudiness at these time scales. One of the 
possibilities would be to model the expected cloud 
amount based on global climate models tuned to 
those past epochs. However, this is not easy to do 
with reliability, as in fact cloud amount and vari- 
ability poses one of the most complicated puzzles 
for climate change today. For example, despite re- 
sults of numerous general circulation models, it is 
unclear how increasing atmospheric temperatures 
resulting from anthropogenic f orcings may influ- 



ence global cloud properties ( Cess et al. 19961) . 



which is important because due to the large ra- 
diative influence exerted by cloud cover, a small 
change in cloud amount or distribution may po- 
tentially provide a strong feedback effect, signifl- 
cantly en hancing or mitigating the effects of globa l 
warmmg (|Dessleijl2010tlSpencer fc BraswelJ201ol ) . 
Here, we use a different approach to tackle this 
problem, by using a semi empirical model to map 
how clouds behave over Earth's surface depending 
on latitude and surface types - such as ice, water, 
desert, and vegetation - which then we apply to 



past epochs. 

To reach this objective the flrst step has been 
to calculate the 1984-2006 climatology of ISSCP 
cloudiness data. With this climatology, we have 
performed a classiflcation of the amount of clouds 
depending on surface types. That is, cloudiness 
data have been split up into four groups according 
to the surface type over which the cloud is located. 
To carry out this classiflcation we have made use 
of real geographical information about the differ- 
ent kinds of surfaces and vegetation in our planet. 
The ISCCP classifies the different surfaces as wa- 
ter, rain forest, deciduous forest, evergreen forest, 
grassland, tundra, shrub-land, desert, and ice, but 
for our study we have only distinguished between 
water, desert, ice, and vegetation, where the latest 
is defined as such zones where there were rain for- 
est, deciduous forest, evergreen forest, grassland, 
tundra, or shrub-land (separate analysis for these 
subgroups gave very similar results). 

Once we got the classification of the cloudiness 
into these four groups, we represented the fraction 
of clouds as a function of latitude for each surface 
type. We have used this scheme to calculate the 
mean cloudiness at each latitude point and its cor- 
responding standard deviation, thus obtaining em- 
pirical relationships between the amount of clouds, 
surface type and latitude. As can be seen in Fig- 
ure [21 these derived cloudiness functions have a 
particular shape which is different for each sur- 
face type, suggesting that global cloudiness dis- 
tribution can be empirically traced according to 
the latitude and the underlying surface type. We 
have repeated this analysis for separate years and 
in seasonal and monthly climatologies with very 
similar results, most of them not shown in this 
paper for space reasons. 

In the next step, these relationships are used in 
an attempt to reconstruct the global mean cloudi- 
ness of the present Earth and the possible cloud 
distribution of past epochs of Earth's history. In 
order to do that, we have taken each surface map, 
both for present and past epochs of the Earth 
(those calculated in Section 2), and we have as- 
signed to each gridpoint a cloudiness fraction by 
using the information of the previously derived 
cloudiness relationships, taking into account both 
surface type and latitude. In this way we obtained 
global mean reconstruction of cloudiness for each 
surface map. 
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As can be seen in Figure [2l because of the con- 
tinental distribution of our own planet, we only 
have information about the fraction of clouds lo- 
cated over vegetation in the latitude range [-60°, 
80°]. Thus in order to reconstruct the cloudiness 
of historical epochs with vegetated areas out of 
this latitude range, we would have to extrapolate 
our functions to get the corresponding cloudiness. 
Since data obtained from extrapolation are sub- 
ject to greater uncertainty, we decided to take the 
cloudiness value corresponding to the higher avail- 
able latitude of the Northern Hemisphere and as- 
signed it as the fraction of clouds of the remaining 
Northern latitudes for which we had no informa- 
tion. For the Southern Hemisphere latitudes for 
which we had no information about cloudiness, 
we assigned to these the fraction of clouds cor- 
responding to the same latitudes in the Northern 
Hemisphere. This filling process does not affect 
substantially our results because the polar regions 
have a low contribution to the total albedo. 

A similar problem arises with clouds located 
over deserts (see Figure [2]). As deserts on the 
present day Earth are mainly located at low lat- 
itudes, in order to avoid extrapolations we have 
chosen to use only two cloudiness values. That 
is, we calculated the mean cloudiness over deserts 
in the latitude range [0° - 30°] and in the range 
[30° - 90°] and we assigned these mean values to 
those points which were located in the respective 
latitude range. These values are 33% for the for- 
mer range and 55% for the last one. Note that we 
have only used this modification in our empirical 
cloud reconstruction for the Cambrian (500 Ma 
ago), since for the other epochs desert areas are 
located mainly at low latitudes. 

In the top left panel of Figure [3l the real Earth 
cloud distribution is shown. In the middle left 
panel it is shown the global cloud distribution 
as reconstructed from our simple model. As can 
be seen in the figure, our semi-empirical model 
for clouds reproduces well the general features of 
Earth's cloudiness distribution. The latitudinal 
range 40°-75°, both north and south, is quite well 
reproduced, as well as the differences in cloud 
amount between land and water at the same lati- 
tude. There are, however, some small-scale struc- 
tures, mostly over tropical oceans, that our model 
does not reproduce, being the major differences 
located in the Pacific and in the Indian oceans. 



These differences are due to oceanic effects, such as 
surface temperature, trade winds, ocean currents, 
and regional meteorological phenomena such as 
"El niiio" / "La niiia" , which can greatly influence 
the cloud formation. This is difficult to parame- 
terize in a general model valid for any continental 
distribution. 

However, with the aim of improving our result 
over the oceanic areas, we have tried to capture 
the longitudinal variations in cloud amount by in- 
troducing an oceanic longitudinal variation in our 
model. We observe that in all large ocean basins, 
although the effect is larger for the Pacific, clouds 
tend to accumulate over the eastern margin (cold 
currents) and have a minimum over the western 
margins (hot currents). Taking the Pacific as the 
rule, we have selected the region where the major 
differences take place, i.e., the region between -20° 
and 0° latitude and between -70° and -180° lon- 
gitude, approximately. We calculated the mean 
cloudiness over this box as a function of longitude 
and fit it with a two-degree polynomial curve (a 
parabola). We have applied this parabolic effect 
to all tropical ocean basins to make a second or- 
der correction to our modeled cloud distributions. 
To do that, the parabola was normalized to the 
mean cloud amount and then the cloud amount in 
each ocean basin was multiplied by this normal- 
ized parabola. As mentioned previously, the same 
correction is applied to all ocean areas within -20° 
to 0°, although the effect is almost unnoticeable 
for small oceans. 

Figure |3] bottom left shows the reconstructed 
cloud distribution of the Earth today after apply- 
ing the longitudinal correction. As can be seen 
cloudiness at the Pacific area has a similar aspect 
to that shown in the same area in the real cloud 
map. In general, a differential map between the 
real cloud amount and our reconstructed cloudi- 
ness clearly indicates that this second order cor- 
rection greatly improves the prediction errors. 

4. LIGHT CURVE RECONSTRUCTIONS 

Once we have been able to reconstruct the 
cloud distribution for present and past epochs of 
the Earth, we want to transform this informa- 
tion into Bond albedo values to test the photo- 
metric variability of the Earth as seen from a dis- 
tant observer. To this end, we have used a sim- 
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pie Earth reflectance model, which calculates the 
visible reflected light curves of our planet for a 
given day, by using as inputs ground scene mod- 
els from the Earth R adiation Budget Experiment 
( Suttles et al cloud and snow/ice cover 

maps, and surface maps as inputs (Pall et al. 
2003). With this information the model is able 
to determine the Bond albedo at each time of the 
day and in any particular direction. Nevertheless, 
the model can only be considered as a flrst-order 
approximation since it does not take into account 
any other climate parameters beyond snow and ice 
that might contribute to changes in Bond albedo. 
Further, the model contemplates only 12 differ- 
ent surface scenes and 4 cloudiness levels (0 — 5%, 
5 - 50%, 50 - 95%, and 95 - 100%). 

Our model calculates first the albedo of each el- 
ement of area of the planet's surface taking into ac- 
count surface type, cloud amount, snow/ice cover, 
and solar zenith angle. Then the model performs 
the calculation of the Bond albedo in terms of the 
albedo of each element of area by integrating over 
all portions of the globe illuminated by the Sun. 
For more details see ( Palle et al., 2001. 

Obviously both cloudiness and surface type 
play a role in the albedo variability and decon- 
volving the effects of both components can be an 
impossible t ask fo r extrasolar planets. However, 
IPalle et ah ( 2008 ) demonstrated that if sufficient 
signal to noise can be obtained, it is possible to 
detect both the presence of continents and the 
clouds. In any case, our aim here is to character- 
ize the photometric variability of both components 
combined. 

4.1. Present Day Earth 

As representative light curves, we have arbitrar- 
ily chosen the months of January, March, and July 
2000 to carry out the Bond albedo simulations of 
the Earth at present day, and at 90, 230, 340, and 
500 Ma ago. In the right panels of Figure |3l we 
plot the modeled 24 hr Bond albedo variation for 
the whole Earth over one day. 

As can be seen in Figure [31 the agreement be- 
tween the light curve obtained from reconstructed 
clouds (bottom right) and that calculated from 
real cloud cover (top right) is quite good. The 
general shape of the light curves are quite simi- 
lar, with the local maximums/minimums located 



at the same hour of the day in both cases. But the 
result obtained from the reconstructed cloudiness 
tend to be higher than that obtained from real 
cloud cover (0.325 as opposed to 0.315). Further- 
more, the amplitude of the light curve is similar 
to that obtained from real cloud data (4.16% and 
3.29% for reconstructed and real cloudiness, re- 
spectively). These differences in the variability of 
the Bond albedo are owing to the fact that real 
cloud cover presents larger changes in the longitu- 
dinal direction in contrast to the low variability in 
such direction of the reconstructed clouds, affect- 
ing Bond albedo values. 

In Figure [3l it can also be seen that the re- 
sulting light curve after applying the longitudi- 
nal correction, i.e., after applying the aforemen- 
tioned parabolic smoothing technique applied over 
oceans, in the Pacific ocean (bottom-right) is more 
similar to that obtained with real cloud data, than 
the light curve calculated before the correction 
(middle- right), confirming that the Earth cloudi- 
ness behavior is better reproduced by using this 
longitudinal correction. 

4.2. Historical Epochs 

In the right panels of Figure lU the light curves 
for the different epochs (90, 230, 340, and 500 Ma 
ago) are shown. As can be seen, the light curves 
of the first three epochs are quite similar. This 
is related to the fact that this epochs have a sim- 
ilar ocean-land-vegetation distribution, having a 
nearly totally clustered continents located princi- 
pally in one hemisphere and deserts located at low 
latitudes. These light curves have a soft shape 
with an enhancement in brightness produced by 
continental masses. This smooth behavior is also 
due to the longitudinal relative homogeneity of the 
reconstructed cloud maps. 

The reflected fight of the Earth at 500 Ma ago, 
however, presents much more variability in con- 
trast to the light curves obtained for the other 
epochs. That is again related to the continental 
distribution. In this epoch, most of the continents 
were clustered in one hemisphere, but it also has 
several big islands that cause strong variability in 
the reflected light. Moreover, this epoch presents 
a significantly higher mean albedo value. That 
can be related to the fact that in this epoch con- 
tinents were covered by deserts, thus involving a 
higher reflectivity. 
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The mean Bond albedo values and their vari- 
ability obtained for the Earth at different epochs 
are summarized in Table [1] By comparing the re- 
sults shown in Table [T] one can note that daily 
variations were larger than present day's, for in- 
stance 500 Ma ago, when the variations were larger 
by a factor of three. Not only that, but also the 
mean albedo values were larger, which should have 
profound influences on the global climate by intro- 
ducing a cooling effect, probably compensated by 
the increased greenhouse gas concentrations over 
that perio d, 16 times that of p re-industrial mod- 
ern levels (jAugustin et al. 



5. CONCLUSIONS 

In this paper, we have used information about 
the surface properties and continental distribution 
of the Earth with the aim of studying the behav- 
ior of the large-scale cloud patterns. We have ob- 
tained empirical relationships between the amount 
of cloudiness and their location on Earth's sur- 
face depending on latitude and underlying surface 
types. This relationships have been used to recon- 
struct the possible cloud distribution for different 
epochs of the Earth history such as 90, 230, 340, 
and 500 Ma ago. With this information, and with 
the help of an albedo model we have attempted, 
for the first time, to reconstruct and understand 
the photometric variability of these past epochs 
according to their different geographic and cloudi- 
ness distribution. 

We find that our model reproduces well the ma- 
jor features of the cloud distribution and the pho- 
tometric light curve of the Earth at present. When 
applied to past epochs of the Earth, we find that 
both the mean albedo value and the diurnal light 
curve variability remain stable as long as desert 
area are confined to the tropical regions. When 
this condition is not met, as during the Late Cam- 
brian about 500 Ma ago, both the mean albedo 
value and the photometric variability are greatly 
increased. This increased variability could help 
in the determination of the rotational period of 
the planet from an astronomical distance. Due to 
the large compositional and chemical changes of 
Earth's atmosphere, we have not attempted to re- 
construct cloud cover maps for epochs prior to the 
Late Cambrian. However, it is likely that the con- 
ditions for this period, i.e., higher albedo values 



and photometric variability, hold for much of the 
previous epochs of the Earth, not considering pos- 
sible albedo variations due to atmospheric changes 
coming from clouds, aerosols, and hazes. 

Research by E. Palle is supported by a Ramon 
y Cajal fellowship from the MICIIN. 
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Table 1: Percentage of the Daily Variability and 
Mean Albedo of the Light Curves for Each Epoch. 



Epoch 


% 


(Albedo) 


Present Earth (real clouds) 


3.29 


0.315 


Present Earth (reconstructed clouds) 


4.16 


0.325 


Late Cretaceous (90 Ma ago) 


4.27 


0.331 


Late Triassic (230 Ma ago) 


5.02 


0.327 


Mississippian (340 Ma ago) 


4.46 


0.329 


Late Cambrian (500 Ma ago) 


12.2 


0.351 



Note. — These quantities have been calculated by performing the mean of these results obtained for 2000 January, March 
and July. 
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Fig. 1. — Global views of the Earth's continental distribution during the Late Cretaceous (90 Ma ago; top 
left), the Late Triassic (230 Ma ago; top right), the Mississippian (340 Ma ago; bottom left), and the Late 
Cambrian (500 Ma ago; bottom right). Courtesy: Ron Blakey, Colorado Plateau Geosystems Inc. 
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Fig. 2. — Empirical relationships between the amount of clouds, surface type, and latitude. These figures 
show the ISCCP 1984-2006 average climatology of cloudiness for each grid cell on Earth as a function of 
latitude (dots), for each surface type separately: water (top left), desert (top right), ice (bottom left), and 
vegetation (bottom right). Solid lines represent the mean cloudiness at each latitude and the error bars 
represent the standard deviation of the mean. 
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Fig. 3. — Cloud distribution and light curves of the present Earth. Left panels illustrate the global distribu- 
tion of real (top) and reconstructed (middle and bottom) cloud cover. Colors represent a different fraction 
of cloudiness as the color code below each panel show. The difference between middle left and bottom left is 
that the last one is the same as the former but after applying the parabolic correction over oceans (see the 
main text). Right panels show the corresponding simulations of Bond Albedo 24 hr variability. The different 
lines styles represent the results obtained for three different months as indicated in the legend. 
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Fig. 4. — Same as Figur^but, from top to bottom, during the Late Cretaceous (90 Ma ago), the Late 
Triassic (230 Ma ago), the Mississippian (340 Ma ago), and the Late Cambrian (500 Ma ago). 
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